Fractional dynamics in the Levy quantum kicked rotor 
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We investigate the quantum kicked rotor in resonance subjected to momentum measurements 
with a Levy waiting time distribution. We find that the system has a sub-ballistic behavior. We 
obtain an analytical expression for the exponent of the power law of the variance as a function of 
the characteristic parameter of the Levy distribution and connect this anomalous diffusion with a 
fractional dynamics. 
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I. INTRODUCTION 

During the last decades it has been possible to ob- 
tain samples of atoms at temperatures in the nano-kelvin 
range [l| (optical molasses) using resonant or quasi- 
resonant exchanges of momentum and energy between 
atoms and laser light. This spectacular experimental 
progress has been accompanied with the development of 
the inter-disciplinary fields of quantum computation and 
quantum information. In this frame, the study of sim- 
ple quantum systems, such as the quantum kicked rotor 
(QKR) 2] and the quantum walk (QW) [1] are extremely 
useful as models to design future codes and computers. 
The behavior of the QKR has two characteristic modali- 
ties: dynamical localization (DL) and ballistic spreading 
of the variance in resonance. These different behaviors 
depend on whether the period of the kick is a rational or 
irrational multiple of 2tt. For rational multiples the be- 
havior of the system is resonant and the average energy 
grows ballistically and for irrational multiples the aver- 
age energy of the system grows, for a short time, in a 
diffusive manner and afterwards DL appears. Quantum 
resonance is a constructive interference phenomena and 
DL is a destructive one. The DL and the ballistic behav- 
ior have already been observed experimentally [1, Q . 

In ref. ^ we investigated the QKR in resonant regime 
and the usual QW when both systems were subjected to 
decoherence with a Levy waiting time distribution. In 
the case of the QKR the model had two strength pa- 
rameters whose action alternated in a such way that the 
time interval between them followed a power law dis- 
tribution. In the case of QW the model used two evolu- 
tion operators whose alternation followed the same power 
law distribution. We showed that this noise in the sec- 
ondary resonances of the QKR and in the usual QW pro- 
duced a change from ballistic to sub-ballistic behavior. 
This change of behavior is similar to that obtained for 
both systems when they are subjected to an aperiodic 



Fibonacci excitation 0, HI ■ In all the above cases the sub- 
ballistic behavior is characterized by the time dependence 
of the variance, i.e. a'^{t) ^ t^", with 0.5 < c < 1. In a 
more recent paper [§| we have studied the QW subjected 
to measurements with a Levy waiting time distribution 
and we found that the system had a sub-ballistic behav- 
ior. We also obtained an analytical expression for the 
exponent of the power law of the variance as a function 
of the characteristic parameter of the Levy distribution. 

In this paper we present a simple model that allows 
an analytical treatment to understand the sub-ballistic 
behavior previously reported in ref. Q. We shall show 
that the temporal sequence of the decoherence, and not 
its intensity, is the main cause of this unexpected dy- 
namics. With this aim we investigate the QKR when 
measurements are performed on the system with waiting 
times between them following a Levy power-law distri- 
bution. We show that this type of noise indeed produces 
sub-ballistic behavior. We obtain analytically a relation 
between the exponent of the variance and the character- 
istic parameter of Levy distribution. These results are 
identical to the ones obtained in ref. Q, showing again 
another aspect of the similarity between QKR and QW, 
as pointed out in previous papers [1, 0, [13, Uli • In ad- 
dition the toy model developed in this work shows that 
a quantum system in combination with a Levy stochas- 
tic process may produce a fractional dynamics for the 
averaged behavior. 



II. LEVY QUANTUM KICKED ROTOR 

The QKR is one of the most simple and best investi- 
gated model whose classical counterpart displays chaos. 
It has the following Hamiltonian 
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where P is the angular momentum operator, / is the 
moment of inertia, K is the strength parameter, 9 is 
the angular position. The external kicks occur at times 
t = n with n integer and unity period. In the angu- 
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lar momentum representation, P\t) = the wave- 

vector is {^{t)) — YlT=-oo^t[tW) and the average en- 
ergy is E{t) = mH\^) = eYZ-oo^^M)\\ where 
e — 121. Using the Schrodinger equation the quantum 
map is readily obtained from the Hamiltonian ((T|) 



a£(t+ 1) 
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where the matrix element of the time step evolution op- 
erator U{k) is 



(3) 



Jjn is the TOth order cylindrical Bessel function and its 
argument is the dimensionless kick strength k = K/h. 
The resonance condition does not depend on k and takes 
place when the frequency of the driving force is commen- 
surable with the frequencies of the free rotor. Inspection 
of Eq. ^ shows that the resonant values of the scale pa- 
rameter T = e/h are the set of the rational multiples of 
27r, i.e. t = 27r p/q. When p/q is an integer the res- 
onance is called principal and when it is a non integer 
rational it is called secondary. 

The dynamics of the Levy quantum kicked rotor 
(LQKR) will be generated by a large sequence of two 
time-step unitary operators Uq and Ui as was done in a 
previous work 6]. But now Uq is the "free" evolution 
of the QKR in resonance and Ui is the operator that 
measures the angular momentum of the QKR. The time 
interval between two applications of the operator Ui is 
generated by a waiting-time distribution p{T), where T 
is a dimensionless integer time step, see Fig. [TJ The de- 
tailed mechanism to obtain the evolution is given in Q. 
We take p{T) in accordance with the Levy distribution 
that includes a parameter a, with < a < 2. When 
a < 2 the second moment of p is infinite, when a = 2 the 
Fourier transform of p is the Gaussian distribution and 
the second moment is finite. Then, this distribution has 
no characteristic size for the temporal jump, except in the 
Gaussian case. The absence of scale makes the Levy ran- 
dom walks scale-invariant fractals. This means that any 
classical trajectory has many scales but none in particu- 
lar dominates the process. This distribution appears, for 
example, in quantum optics [Tsj as an appropriate tool 
to describe cooled atomic samples in terms of a compe- 
tition between a trapping process (the atom falls in the 
optical trap) and a recycling process (the atom leaves the 
trap and eventually return to it). The most important 
characteristic of the Levy noise is the power-law shape of 
the tail, accordingly in this work we use the waiting-time 
distribution 
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< t < 1, 

t > 1. 



(4) 



To obtain the time interval T we sort a continuous vari- 
able t in agreement with Eq. ^ and then we take the 
integer part Ti of this variable [Q] . 
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FIG. 1: Paths of the LQKR wave function as a function of 
the angular momentum, I. The measurements are performed 
at times Ti and the wave function collapses at these times. 
Between measurements the system has an unitary quantum 
evolution 



In what follows we assume that the resonance condition 
of the QKR is satisfied, for the sake of simplicity we take 
e/h — 27r in such a way that the operator Uq corresponds 
to the first principal resonance. This choice does not 
imply a loss of generality for our results as we shall show 
below. 

Let us suppose that the wave function is measured at 
the time t, then it evolves according to the unitary map 
Eq. ([2]) during a time interval T, and again at this last 
time t + T Si new measurement is performed. In Fig. [1] we 
present a path diagram of the state evolution. It shows 
four time steps when the measurements are perform, be- 
tween measurements there is an unitary evolution. When 
the measurement is performed the wave function col- 
lapses in a momentum state. The resulting states after 
successive measurements need not be contiguous states 
as in the QW because all transitions are possible. 

In the figure we present a generic and arbitrary path 
with bold line. From this diagram we can write a dynam- 
ical equation for the probabilities of the LQKR momenta. 
To begin note that the probability that the wave function 
collapses in the eigenstate due to a momentum mea- 
surement, starting from the eigenstate |0) after a time T 
is 



P,{T)^\a,{T)\ 



(5) 



The momentum distribution Pj depends on the initial 
state and on the time interval T because of the collapse 
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of the wave function, and it will play the role of transi- 
tion probabilities for the global evolution. The mecha- 
nism used to perform momentum measurements assures 
that these distributions will repeat themselves around the 
new momentum. Then it is straightforward to build the 
probability distribution Pj at the new time t + T as a 
convolution between this distribution at the time t and 
the conditional probability: 



m+T)= J2 qi-,{T)p,{t), 



(6) 



] = -oc 



where qi^j are the transition probabilities from state j to 
state I and the sum is extended between — oo and -l-oo be- 
cause all the transition are possible. To calculate qi^j the 
original dynamical equations Eqs. ^ [3]) and the proper- 
ties of the Bessel function are used to obtain a connection 
between the initial pure state after a measurement and 
all possible final states before the next measurement 



(7) 



The Eq. ([6]) is a sort of master equation, but not strictly 
because of the time dependence of the transition proba- 
bilities. 

There are many ways of solving Eq. ^ , we choose the 
method of the generating function G(z.t) defined as 



G{z,t)= £ z^P,{t), 



(8) 



j = -oo 



where we shall take the auxiliary variable z as z = e^^ , 
with if real. It is easy to prove using Eq. ([5]), that 

G(z, t + T) = G{z, t)jQ{2KTsinip/2), (9) 

The first and second moments are given by 



mi(t) =G'(1,<), 
m2(i) = G"(1,<), 



(10) 
(11) 



where the prime indicates differentiation with respect to 
z. Then using these equations and Eq. ^ the following 
maps for the moments are obtained: 

mi(i + T) = TOi(t)+mi,(T), (12) 
m2{t + T) ^ m^it) + 2mi{t)miq{T) + m2g{T), (13) 



where 
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(14) 

m2,(r) = }^ I'qiiT). (15) 

— oo 



Note that miq{T) and m2q{T) are the first and sec- 
ond moments of the unitary evolution between measure- 
ments. From these expressions and using the Eq.([7]) 



the following results are obtained miq(T) = and 
m2q{T) = j^r^. Therefore the global variance cr^(t) = 
m2 it) — m\ (t) verifies that 



(16) 



where a^iT) = m2q{T) - ml^(T) = \k,T'^ is the vari- 
ance associated to the unitary evolution between mea- 
surements. Note that the value of the coefficient of 
is a consequence of using the principal resonance but the 
time dependence remains unchanged for any other reso- 
nance. From these last equations is easy to show that 



1 ^ 



where 



N 



(17) 



(18) 



and N is the number of measurements performed. These 
results are generic, now we shall calculate the average of 
Eq. (HZl) 



1 (T^) 
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where the relation t ~ {Ti){N) was used. The first and 
the second moments of the waiting time for our Levy 
distribution, Eq. (|4]), are 



^|(H^),"^l- (20) 
a + l\ (i+ln(t)), a=l. 



1 I t^~°-i 

3 ^ 2- 



a+1 [ (|+ln(t)), a = 2. 



(21) 



Substituting these expressions in Eq. p9|) and for a large 
time 



1 



('^'W)=2^r(3-a)^ if 1 



, if a 1, 



(22) 



Therefore when t ^ oo the variance behaves as (cr^ (t)) 
t^'^ where 



(3 - a) 



if < a =^ 1, 
if 1 a 2. 



(23) 



This result shows that measurements do not break com- 
pletely the coherence of the system on a time scale that 
includes several of them. For ^ a ^ 1 the ballistic 
behavior is preserved as in the usual resonant QKR, and 
for 1 ^ a ^ 2 it is lost and the sub-ballistic behavior 
takes place. When a = 2 the system has a diffusive be- 
havior as in the usual Brownian motion. From the fact, 
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that the exponent c does not depend on k it follows that 
the results are valid for both primary and secondary res- 
onances. However the coefficient c depends on the micro- 
scopic law of evolution, o-giT) oc T^. Then, we may pose 
the question if there exists a relation between the time 
dependence of cr^(T) in the quantum unitary evolution 
between measurements and the exponent c of the power 
law for the averaged variance. To answer this question we 
shall suppose a unitary quantum evolution that produces 
the following variance 

al{T)^T^, (24) 

where /? is a constant. The reasoning to obtain the ex- 
ponent c can now be repeated, it is only necessary to 
calculate again the new expression for a general (3 mo- 
ment with the Levy waiting time distribution, that is 

\ / « + (^+ln(t)). a = (3 

Then, in this generic case, for t oo, the exponent c is 

{(3, if a < /3,0 s$ a 1, 

(/3-a + l), if «</3,l s;as;2, 

a, if /3 =s; a,0 s; a 1, 

1. if/3^a,l<Q!<2 

This expression shows that these systems can exhibit dif- 
fusive, sub-diffusive, ballistic or sub-ballistic behaviors 
depending on the values of a and /3. Then, in the theo- 
retical frame of fractional dynamics 14] Eq. ([6]) together 
with the Levy distribution would generate a generalized 
master equation from which a generic fractional diffusion 
equation 'T§\ could be built. This fractional dynamics 
approach has as an extreme case the classical diffusion 
equation for a — /3 — 2. 



with a far-detuned standing wave of laser light [j, Il6l |. 
We study the QKR subjected to measurements with a 
Levy waiting-time distribution. As the Gaussian distri- 
bution is a particular case of the Levy distribution, then 
our study is open to wider experimental situations. We 
showed numerically [6] that a Levy noise does not break 
completely the coherence in the dynamics of the QKR 
but produces a sub-ballistic behavior. There the system 
was also a LQKR but the operators Uq and Ui corre- 
sponded to the same secondary resonance with two dif- 
ferent values of the strength parameter k and these oper- 
ators do not commute. It is important to note that if the 
operators correspond to a primary resonance the ballistic 
behavior was retained due to the commutativity between 
the operators Uq and Ui In the present model one of 
the operators is unitary, and may correspond to any reso- 
nance of the QKR, and the other is the measurement op- 
erator. These operators also do not commute and again 
this is linked to the sub-ballistic behavior. Then we can 
conclude that for the LQKR the behavior of a(t) depends 
on the commutativity and the waiting time distribution, 
both models show the same physics. We developed a 
simple analytical theory to connect the waiting time pa- 
rameter a with the exponent c. The LQKR behaves like 
the QW subjected to the same measurement process 
strengthening our previously established parallelism be- 
tween both systems [1, 0, IHj, where the resonant 
QKR is interpreted as a QW in momentum space. The 
type of model developed in this work shows that a quan- 
tum system in combination with a Levy stochastic pro- 
cess leads to an anomalous diffusion and not to the well 
known diffusive process of Browniam motion. Finally, 
this simple toy model may help to understand the con- 
nection between a fractional approach and a generalized 
master equation. 



III. CONCLUSION 

The quantum resonances of the QKR have been exper- 
imentally observed in samples of cold atoms interacting 
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